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,he  aUlhor5  formulate  .he  optimum  replacement  problem  as  a  concave 
integer  programming  problem,  and  develop  an  algorithm  to  solve  the 
formulated  problem. 
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Formulation  of  the  Problem 


'.v'e  shall  consider  a  system  consisting  of  n  elements.  Assume  that  the 
probability  of  failure  in  a  definite  time  period  is  known  for  each  element 
and  equals  q^.  In  addition  to  the  basic  element,  there  are  Zj  spare 
elements  (z^  is  the  number  of  replacements) ,  It  is  assumed  that  the  spare 
elements  function  in  parallel  with  the  basic  element,  or  that  they  can  be 
instantaneously  connected  with  the  system  when  the  main  element  fails. 

A  particular  element  of  the  system  is  considered  to  have  failed  during 
the  one-time  failure  of  the  basic  element  and  all  spare  elements.  Then  the 
probability  of  failure  during  a  definite  time  period  equals,  for  the  ith 
element,  qi*+zi  where  the  number  of  spares  is  z^. 

Suppose  that  the  failure  of  any  of  the  n  elements  leads  to  the  failure 
of  the  complete  system.  Then  the  probability  that  the  system  will  not  fail 
during  a  definite  time  period  equals 


(0  «  n  0  -  <*>> 


(1) 


where  »  1  ♦  z±. 

The  quantity  ♦  can  be  considered  as  a  quantitative  indicator  of  the 
degree  of  reliability  of  the  system.  From  here  on,  for  brevity,  we  shall  call 
this  quantity  the  reliability  of  thi.  system. 

Every  ith  element  of  the  system  is  characterized  by  definite  dimensions: 
weight,  cost,  etc.,  which  are  expressed  in  terms  of  the  coefficient  a^j  (jth 
characteristic  of  the  ith  element).  Assume  that  admissible  dimensions,  weight 


and  other  similar  characteristics  of  the  system  bj ,  are  given.  Then  the 
optimal  replacement  problem  can  be  formulated  as  follows.  To  determine 

max  <D  ■■■-■  n«- 

<-i  4 


with  the  conditions 


n 

S2i  a0xi  <  bh 


my, 


(3) 


where  >  1  are  integers.  (4) 

In  the  context  of  the  problem,  Xj  must  be  integers,  since  these 
represent  the  number  of  elements .  Here  the  optimality  criterion  $  is  a 
concave  function  of  the  variables  xi.  Thus  the  optimal  replacement  program 
is  formulated  as  a  concave  integer  programming  problem.  At  the  present  time, 
direct  methods  for  solving  problems  of  this  class  do  not  exist. 

Problem  (2)  through  (4)  can  be  replaced  by  the  problem  of  finding 
n  n 

min  [-In  $>)  =  min  £  -  ln(l  -  q*i)  -  min  I  <#>i  (x.) 

i=l  1  i-1 

subject  to  restrictions  (3)  and  (4) . 

Here  (x^)  =  -ln(l  -  is  a  convex  function  (see  the  figure,  curve  1) . 

The  objective  function  for  this  problem  has  the  property  that  it  is  a  sum  of 
convex  functions,  each  of  which  depends  only  on  a  single  variable.  A  procedure 
for  reducing  this  problem  to  a  linear  programming  problem  is  described  in  the 
literature.^  The  optimal  solution  for  the  latter  problem  is  an  integer  solution. 


!  D.B.  Yudin  and  Ye.G.  Goldstein,  "Problems  aid  Methods  of  Linear 
Programming,"  Sovetskoye  Radio,  Moscow,  1961.  ! 
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However,  a  practical  application  of  this  method  is  difficult  because  of 
the  high  dimensionality  of  the  corresponding  linear  programming  problem.  The 
same  remark  also  applies  to  dynamic  programming  methods,  the  application  of 
which  to  problems  (2)  through  (4)  has  been  studied.!  •' 

An  important  special  case  of  the  optimal  replacement  problem  is  the  case 
when  the  dimensions  of  the  system  are  restricted  only  to  one  parameter,  for 
example,  weight  or  volume.  Then  the  number  of  restrictions  of  type  (3)  equals 
one  (m  =  1) . 

In  this  case,  the  problem  can  be  solved  efficiently  by  tho  following 
method. 

We  will  consider  the  problem: 

tt 

min  min  FX. 

with 

7i 

li  a\xi  ^  b, 

i-t 


where  X£  >  x^q  can  assume  only  discrete  values 

*io -h  kh,  (*=-1,2,...). 


Here  <J>^(xj)  are  convex,  monotonically  decreasing  functions  of  a  single 
variable;  h  >  0  is  the  discrete  step  size,  b  >  0  and  aj  >  0  (i  «  1,...,  n). 

For  the  all-integer  case,  xig  =  1,  and  h  =  1.  Generally  speaking,  the 
discrete  step  size  can  be  different  for  various  variables:  h  *»  h(i). 

This  case  is  easily  reduced  to  the  case  h  =  const  by  introducing  new 
variables  of  the  type  x'j  »  kx^  and  is  the  case  considered  here. 

Suppose  that  for  all  i,  »  xjq.  It  is  clear  that  in  this  case 
inequality  (6)  is  satisfied,  since  in  the  contrary  case  a  solution  does  not 
exist.  We  will  determine  a  decrease  in  the  objective  function  (5)  by  increasing 
each  of  the  unknowns  by  h:  4>i(xi0)  -  $i(xio  +  h)  >  0(1  =  1,  2,  n) .  We  remember 

that  are  monotonically  decreasing  functions. 

Dividing  this  quantity  by  a^,  we  obtain  6^  --  the  decrease  in  the 
objective  function  per  unit  weight.  In  the  general  case,  by  weight  we  mean  a 


(5) 


(6) 


(7) 


1965. 


!  R.  Bellman  and  S.  Dreyfus,  "Applied  Dynamic  Programming,"  Nauka,  Moscow, 


property  which  can  be  expressed  by  the  coefficient  b  where  the  unit  weight 
of  each  variable  equals  a^. 

Among  all  unknowns  we  find  that  unknown  which  insures  a  maximum  decrease 
in  the  objective  function  per  unit  weight,  i.e.,  we  find 

run  x  — *  ■»  — ———————  bs  «•<- — ■  "  ■  ■  • 

l  n{  a  k 

We  .increase  the  value  of  xj<  by  h  and  verify  condition  (6) .  If  it  is 
satisfied,  then  for  we  have  a  new  value  which  has  been  increased  by  h 
in  comparison  with  the  previous  value.  We  again  compute  the  quantities  (8), 
but  here,  instead  of  xjcg  we  take  the  new  value  x^q  +  h.  We  notice  that 
the  quantity  (8)  must  be  computed  only  for  i  =  k,  because  for  all  remaining 
i  values  it  remains  unchanged.  This  considerably  shortens  the  computational 
effort.  We  again  determine 

_ <r4  (*-,)  -  ¥<  (*,  H-  A)  <pr  (*,.)  -  <pr  (rr  h) 

- — —  ""  =  — — - — —  — . — - 

i  at 

and  if  condition  (6)  is  satisfied,  we  increase  xr  by  h.  The  computation 
continues  until  either  the  inciease  of  the  current  xr  by  h  leads  to: 
a)  attainment  of  equality  in  (6)  (then  the  value  xr  is  replaced  by  xr  +  h, 
and  the  computation  terminates) ,  b)  the  violation  of  condition  (6) . 

We  then  determine  xj ,  for  which  the  quantity 

Vj  (*•;)  —  <Pj  (’•)-;■/»)  _ 

ni  '  1 

assumes  a  value  which  is  next  to 

max 

i 


If,  while  increasing  this  Xj  by  h  condition  (6)  is  satisfied,  Xj 
is  replaced  by  xj  +  h  and  the  calculation  continues.  If  among  all  xj  there 
is  no  xi  which  can  be  increased  by  h  without  violating  condition  (6),  compu¬ 
tation  terminates . 

It  can  be  shown  that  the  method  shown  above  gives  an  exact  solution  to 
problems  (5)  through  (7)  if  it  terminates  in  attainment  of  equality  in  (6) . 

If  for  the  solution  X*  which  was  obtained  condition  (6)  has  the  form 
of  an  inequality,  then  the  solution,  generally  speaking,  need  not  be  optimal. 

In  this  case  let  X0pt  be  the  optimum  solution.  Then  F(X*)  >  F(X0pt). 

We  will  determine  a  solution  x  which  will  be  obtained  when  condition  (6) 
is  violated  for  the  first  time,  and  without  regard  for  this,  the  corresponding  x 


are  increased  by  h  and  the  computation  is  considered  terminated.  Clearly 
then  the  inequality  F(X*)  >  F(X0pt)  * 

From  here  we  have  an  estimate  for  the  constructed  solution 
F(X*)  -  F(Xopt)  <  F(X*)  -  F(X). 

For  the  optimum  replacement,  problem, 


a. 


<Mr<)~<Vr<  i*  h) 


■  la  (1  —  -f  In  (1  —  q***')  _ i_ .  1  -  ' 


In  addition  to  problems  (S)  through  (7) ,  the  above  method  is  applicable  in  the 
following  cases. 


1.  6i(xig)  >  0  and  are  convex  functions  which  are  not 

necessarily  monotonically  decreasing  (curve  2).  In  this  case,  every  time 
after  increasing  the  current  by  h,  one  must  verify  the  sign  of 

«i Cxi  +  h).  In  addition  to  the  usual  condition  for  terminating  the 
calculation,  in  this  case,  each  can  be  increased  only  until 

the  condition  6*  >  0  takes  place. 

2.  <t>i(Xi)  are  monotonically  increasing  convex  functions,  with 
Si'xio)  <  0  (line  3)  and  condition  (6)  has  the  form 


n 


^a4x,>6. 


(6') 


3.  $i(xi)  are  convex  functions  and  condition  (6')  holds.  In  this  case, 
for  XjQ  we  can  take  some  quantities  which  exceed  considerably  their  optimum 
values  and  for  which  6i(xio)  <  0,  and  condition  (6')  is  satisfied.  The  move¬ 
ment  takes  place  in  the  direction  of  decreasing  Xj. 

4.  Problems  (5)  through  (7),  in  the  case  of  continuous  change  of  variables. 


For  a  sufficiently  small  step,  an  optimum  solution  can  be  attained  with 
arbitrary,  preassigned  precision.  To  reduce  the  amount  of  computation,  the 
step  size  need  only  be  decreased  sufficiently  near  the  optimum. 
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